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MORPHOLOGIES OF HIGH REDSHIFT, DUST OBSCURED GALAXIES FROM KECK LASER GUIDE STAR 

ADAPTIVE OPTICS 
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ABSTRACT 

Spitzer MIPS images in the Bootes field of the NOAO Deep Wide-Field Survey have revealed a class of extremely 
dust obscured galaxy (DOG) at z ~ 2. The DOGs are defined by very red optical to mid-IR (observed- frame) 
colors, R — [24 fim] > 14 mag, i.e. /^(24 fim)/ f v (R) > 1000. They are Ultra-Luminous Infrared Galaxies with 
^8-iooo/xm > 10 12 — 10 14 L , but typically have very faint optical (rest-frame UV) fluxes. We imaged three DOGs 
with the Keck Laser Guide Star Adaptive Optics (LGSAO) system, obtaining ~ 0.06" resolution in the i\~'-band. One 
system was dominated by a point source, while the other two were clearly resolved. Of the resolved sources, one can 
be modeled as a exponential disk system. The other is consistent with a de Vaucouleurs profile typical of elliptical 
galaxies. The non-parametric measures of their concentration and asymmetry, show the DOGs to be both compact 
and smooth. The AO images rule out double nuclei with separations of greater than 0.1" (< 1 kpc at z = 2), making 
it unlikely that ongoing major mergers (mass ratios of 1/3 and greater) are triggering the high IR luminosities. By 
contrast, high resolution images of z ~ 2 SCUBA sources tend to show multiple components and a higher degree of 
asymmetry. We compare near-IR morphologies of the DOGs with a set of z = 1 luminous infrared galaxies (LIRGs; 
Lir ~ 10 11 L©) imaged with Keck LGSAO by the Center for Adaptive Optics Treasury Survey. The DOGs in 
our sample have significantly smaller effective radii, ~ 1/4 the size of the z = 1 LIRGs, and tend towards higher 
concentrations. The small sizes and high concentrations may help explain the globally obscured rest-frame blue-to-UV 
emission of the DOGs. 

Subject headings: galaxies: high-redshift — galaxies: structure — infrared: galaxies — instrumenta- 
tion: adaptive optics 



1. INTRODUCTION 

The Spitzer Space Telescope (|Werner et al.ll2004f ) has 

i-H . revealed large numbers of luminous and ultra-luminous 
infrared galaxies (LIRGs with Lir ~ 10 11-12 L©, 

^^ \ and ULIRGs with L T r > 10 1 2 L&) in the distan t 
universe (e.g., iBell et all 120051 iLe Floc'h et~afl l2005h . 

^^ ■ Dust heating by star formation and active galactic nu- 
clei (AGN) within the LIRGs and ULIRGs result in 
large mid-infrared fluxes detectable to high redshift 
by the Spitzer M ultiband Imaging Photometer (MIPS; 
Rie ke et aDl2004f ). Because much of the energy in these 

S^ systems has been reprocessed by dust, optical (rest-frame 
UV) stu dies have underes timated their bolometric lumi- 
nosities (IBell et al.l 12005). At the extreme are sources 
that are highly obscured at observed optical wavelengths. 
Recent MIPS 24 /im images in the Bootes field 
of the NOAO Deep Wide-Field Survey (NDWFS; 
iJannuzi fc Devi [l999) have reveale d a class of ex- 
tremely dust ob scured galaxy (DOG: iHouck et aUl2005l : 
IWeedman etHI 12006b Brand et all I2007T ). These sys- 
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terns were selected to have 24 /im flux densities in ex- 
cess of 0.3 mJy, and optical to IR colors redder than 
R — [24] > 14 mag, rou ghly a magni tude redder than 
Arp 220 at any redshift (|Dev et al.ll2008l ). Spectroscopic 
redshift confirmations for 86 DOGs were obtained with a 
combination of optical spectroscopy (Keck DEIMOS and 
LRIS; iDesai et al . 2007), near-IR spectroscopy (Keck 
NIRSPEC: lBrandet al. 2007), an d mid-IR spectroscop; 
(Spitzer IRS; IHouck et all 120051 : IWeedman et~afl 1200 
Higdon et al., in preparation). The redshift distribution 
of these DOGs is Gaussian, with a mean of z — 1.99, a 
sigma of a* = 0.5, and a redshift range of z — 0.8 — 3.2 
(jDev et al.H2008[ ). 

The rest-frame UV to mid-IR spectral energy distribu- 
tions (SEDs) of DOGs range from power laws, rising into 
the mid-IR, to those that show a "bump" at rest wave- 
length 1.6 \±vn (observed in Spitzer IRAC 3-8 \±m bands; 
Desai et al. in preparation, Dey et al. 2008). A possible 
explanation for these two classes is that the power law 
sources are AGN dominated, whereas the SEDs of the 
"bump" sources are primarily powered by star forma- 
tion. For instance, spectroscopic Ha observations of 10 
power law dominated, dust obscure d galaxies in Boote s 
showed that all 10 harbor an AGN (|Brand et al.ll2007f ). 
In reality, most DOGs are likely to contain some combi- 
nation of bo th vig orous star formation and AGN activity 
(jFiore et alJl2Q08h . 

While spectra and SEDs may reveal the primary power 
sources for the DOGs, they do not necessarily reveal the 
triggering mechanism. Are the DOGs undergoing violent 
mergers that drive both central star formation and AGN? 
Alternatively, are the DOGs isolated massive galaxies 
seen at the time of their first collapse? Are there other 



processes contributing such as bar formation or minor 
mergers? Deep, spatially resolved imaging of the DOGs 
should help distinguish between some of these scenarios. 
For instance, recent merging acti vity may be revealed b y 
multiple nuclei or tidal tails (e.g.. I Melbourne et~a l. 2005). 
Smooth compact systems may be mor e readily explaine d 
by a primordial collapse model (e.g., IZirm et alJl2007l ). 
Imaging may also be useful for identifying strong cen- 
tral point sources suggesting AGN activity. In order to 
observe the morphology of the underlying stellar popula- 
tion, rest-frame optical to near-IR imaging is preferable 
to imaging at shorter wavelengths. It is less affected by 
ongoing star formation and dust obscuration compared 
with rest-frame UV imaging. Rest-frame UV, however, 
is valuable for identifying sites of recent star formation 
and unobscured AGN, which may or may not track the 
underlying stellar distribution. 

We have undertaken two efforts to obtain resolved im- 
ages of the DOGs. This paper reports on high spatial 
resolution (0.06") Keck laser guide star (LGS) adaptive 
optics (AO) if'-band imaging (rest-frame optical to near- 
IR) of three DOGs. While this sample is small, it is 
the first and only set of high resolution K-band images 
available for the DOGs. Bussman et al. (in preparation) 
will report on Hubble Space Telescope (HST) imaging 
of an additional 30 DOGs. The Bussman et al. data 
include ACS and WFPC2 /-band imaging which probe 
rest-frame blue to UV, and which will be valuable for de- 
tecting the sites of ongoing star formation. The Bussman 
et al. data also include NICMOS i/-band imaging (rest- 
frame optical), which will be more sensitive to the un- 
derlying stellar population. In comparison to NICMOS, 
the Keck AO data presented in this paper probe longer 
wavelengths, and have roughly a factor of 3 higher spa- 
tial resolution, providing an unprecedented view of the 
stellar distributions in the cores of the DOGs. 

Section 2 describes the sample selection, the observa- 
tions, and the data reduction. Section 3 describes the 
morphologies of the DOGs and provides measurements 
of their structural parameters. A discussion and com- 
parison of the DOGs with two additional samples of high 
redshift dusty galaxies is given in Section 4. Through- 
out, we report Vega magnitudes and assume a A cold- 
dark-matter cosmology: a flat universe with Hq = 70 
km/s/Mpc, Q m = 0.3, and Q\ = 0.7. 

2. THE DATA 

Keck LGS AO observations of three DOGs were ob- 
tained on the night of U.T. 2007 May 21. The AO system 
tracks and corrects wavefront distortions of astronomical 
sources as the light propagates through the Earth's tur- 
bulent atmosphere. For a complete descr i ption of the 
Keck LGS AO system, see Ivan Dam et al.l (|2006[ ). De- 
spite the laser, the AO system still requires a nearby 
(< 55") faint (R < 18) natural guide star ("tip-tilt" 
star) to correct for image motion. By correcting for the 
effects of turbulence, the AO system achieves diffraction 
limited spatial resolution. The Keck AO system has a K- 
band spatial resolution of ~ 0.06", comparable to HST 
in the optical, and a factor of ~ 3 higher resolution than 
HST NICMOS in the near-IR. The performance of the 
AO system can be sum marized by the Strehl ratio (e.g., 
IWizinowich et aLll200QJ ). This is the ratio of the peak 
brightness of the observed point-spread-function (PSF) 
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Fig. 1. — Spectral energy distributions (SEDs) of the three DOGs 
in our sample. Photometric measurements were made from opt i- 
cal NDWFS B W ,R,I ground based images (Jannuzi & Dev 1999), 
near-IR J, Ks FLAMEX images ( Elston et al. 2006) , the four 
Spitzer IRAC channels (3.6, 4.5, 5.8, and 8 /im Eisenhardt et al. 
2004), and the Spitzer mid-IR 24 fim band (Houck et al. 2005). 
DOGs 1 and 2 show power-law SEDs, while the SED of DOG 3 ap- 
pears to have a "bump" in the IRAC bands. Typical photometric 
uncertainties are roughly the size of the points. 
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Fig 2. The redshift distribution of the DOGs fBlack [Dev et al.l 
2008). Also shown is the distribution for different X-band magni- 
tude limits. The K-band bright targets (K < 19, blue, top panel) 
favor lower redshifts (z ~ 1) but still span the range of the full 
redshift distribution (Black). In contrast, fainter K-band sources 
(green, middle panel) trace the larger distribution more closely. A 
significant number of the DOGs with redshifts were either not ob- 
served at K, or were too faint to be detected (red, bottom panel). 

to the theoretical maximum peak of a perfect diffraction 
limited PSF. Table [1] gives a summary of the observa- 
tions, including exposure times and AO performance as 
given by the full-width-half-maximum (FWHM) of the 
PSF and the Strehl ratio. 

2.1. The Sample 

We identifie d a set of 2603 D OGs in the Bootes field of 
the NDWFS (jDev et all 12008). These galaxies were se- 
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Fig. 3. — A color-magnitude diagram of the DOGs found in 
the NDWFS (black points). DOGs with spectroscopic redshifts 
are circled in blue. The AO sample is shown with red diamonds. 
There does not appear to be anything unusual about the DOGs in 
the AO sample. 

lected to have R- [24] > 14, roughly ^(24 /am)/f v (R) > 
1000, and 24 /im flux densities in excess of 0.3 mJy. We 
matched these objects to the positions of suitable tip-tilt 
AO g uide stars from the USNO-A2.0 catalogue (|Monetl 
1998). For this pilot program, we selected DOGs that 
were located within 30" of stars brighter than R = 15.5, 
optimal conditions for achieving the highest level of AO 
correction. This provided us with a set of 55 potential 
targets. 

We prioritized this sample by measured ifs-band 
mag nitude (from the FLAMEX survey of Bootes; 
lElstonet al J 12003 ). For DOGs with no FLAMEX Ks 
imaging, Ks magnitudes were extrapolated from Spitzer 
IRAC 3.6 /am photometry. Due to poor weather, we were 
only able to image three DOGs (denoted DOG 1, DOG 
2, and DOG 3 for simplicity). In order to ensure suc- 
cessful observations of a preliminary sample despite lim- 
ited observing time, the three DOGs selected were among 
the brightest ifs-band sources in the sample (Ks < 19). 
Note: while DOG 3 did not have a FLAMEX observa- 
tion, its extrapolated if s-magnitude was brighter than 
this limit. Subsequent Ks observations of DOG 3, made 
by our tea m with the Palomar Wide-field Infrared Cam- 
era fWIRC: IWiTson et alJl2003h . confirmed this. With 24 
jam fluxes > 1 mJy, they also happened to be among the 
brightest 24 jam sources in our set of 55 possible targets. 
The MIPS 24 /am fluxes and the FLAMEX and WIRC 
K 5-band magnitudes are given in Table O 

A final consideration in the sample selection was given 
to the shape of the optical - IR SED of the DOGs. 
We chose our sample such that it would contain at 
least one of each SED type, power law, and "bump". 
Figure [T] shows the SEDs of the three DOGs in our 



sample. The ground-based optical dat a, B w , R, and 
/-ban d filters, are from the NDWFS (|Jannuzi fc Devi 
1999). The ground-based nea r-IR data, J and if s-band 
filters, are from FLAMEX (jElston et all l2QQ6h . The 
four space-based Spitzer IRA C channels, 3.6, 4 . 5, 5.8 , 
and 8 /am, were described bv lEisenhardt et all J2004), 
while the Sp itzer mid-IR 24 /am data were described by 
iHouck et al.l (|2Q05h . Two of the galaxies selected, DOGs 
1 and 2, show power law SEDs. DOG 3 is the most 
likely "bump" candidate. While we succeeded in select- 
ing at least one of each SED type, we caution that DOG 
3 is not the best "bump" candidate in the larger DOG 
pool. There are other DOGs in Bootes that show a more 
prominent "bump" feature. 

While DOG 2 appears to have a power law SED, the 
if s-band data point is offset from the rest of the power 
law. This data point is based on the FLAMEX K s-band 
image, and is measured to within 10%. In addition, we 
checked the FLAMEX zeropoint with 2MASS images of 
the field and found good agreement. It is possible that 
a strong emission line is contributing to the deviation of 
the K s-band data point from the power law SED. If the 
object is at z ~ 2 then Ha falls into the Ks pass-band. 
A very large, rest-frame Ha equivalent width, on the or- 
der of 1000 A, could account for this offset. This is about 
twi ce as large as t he lar gest Ha equivalent width found in 
the lBrandetafl ([2007h sample of NIRSPEC DOG spec- 
troscopy. While this data point remains puzzling, the 
remainder of the analysis in this paper is independent 
of the total measured if -band flux. Instead we will be 
investigating how the flux is distributed spatially. 

Unfortunately, because of their proximity to bright AO 
guide stars, the three DOGs in the AO sample have not 
previously been targeted by our spectroscopic surveys. 
Whil e we do not know the exact redshifts of these galax- 
ies, iDey et all (|2008l ) has shown that the DOG selection 
criteria strongly selects for z > 1 dust obscured galax- 
ies. Figur e [2] sho ws the redshift distribution as given in 
iDev et all (|2008h . Also shown are the redshift distribu- 
tions at different if -band flux levels. While the if -band 
bright DOGs show a preference for the low-z end of the 
redshift distribution (z ~ 1), they do span the range of 
the larger DOG sample. As we primarily limit our anal- 
ysis to morphology, a lack of redshift information does 
not prove critical to our analysis. 

Figure [3] shows how the photometry of the AO sample 
compares to the larger DOG sample. Aside from being 
at the brighter end of the distribution, there does not 
appear to be anything unusual about the three DOGs in 
the AO sample. 

2.2. AO Observations 

AO observations were made in the if '-band with the 
NIRC2 camera. In order to achieve the full Keck reso- 
lution, we observed in the narrow field mode which has 
an 0.01" pixel scale, and a 10" field of view. Individual 
exposures were 180s and a dither was applied between 
exposures. DOGs 1 and 2 were observed for 30 minutes 
each, while DOG 3 was observed for 9 minutes. Obser- 
vations of DOG 3 were cut short by a laser fault. 

In Section 3 we model the two-dimensional light pro- 
files of our galaxies. In order to do this accurately, we 
tracked the complicated AO PSF throughout the night. 
The 10" x 10" field surrounding DOG 1 contained a faint 



TABLE 1 

Keck Observation Summary 



object MIPS Catalogue Name a exptime FWHM Strehl Guide Star 

[min] of PSF ["] estimate [%] R [mag] sep ["] 



DOG 1 SST24 J143234.9+333637 
DOG 2 SST24 J142801. 0+341525 
DOG 3 SST24 J142944.9+324332 
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Fig. 4. — K' images of the DOGs 1 (top), 2 (middle), and 3 (bottom). The left column gives the observation. The middle column gives 
the best fit model of the light distribution. The right column is the residual difference between the two (observation - model). DOG 1 is 
best fit by an exponential disk with Sersic n = 1.5. DOG 2 is best fit by a de Vaucouleurs profile (actual n = 3.5). DOG 3 appears to be 
dominated by a point-source. The images of the DOGs tend to be smooth showing little evidence for substructure, such as multiple nuclei. 
The faint substructures seen in the residual images, especially for DOGs 1 and 3, represent a small fraction of the total flux. At the center 
of the galaxy, the RMS deviation of the residuals is 5% or less of the galaxy flux/pix. The images are ~ 0.6" on a side, with north up and 
east to the left. 



star suitable for tracking the real time variations of the 
AO PSF. We used this PSF to model the spatial light 
profile of DOG 1. In order to track the PSF of DOG 2 we 
obtained images of the tip-tilt star immediately follow- 
ing the science observations. Because of the laser fault, 
we did not obtain a PSF specifically matched to DOG 3, 
but use both existing PSFs in the analysis of the DOG 
3 light profile. 

2.3. Data Reduction 

The observations were reduced in a manner similar 
to iMelbourne et al.l (|2QQ7l ). Frames were corrected for 
flatfield variations with an average combined dome flat. 
Frames were then sky-subtracted using a scaled clipped 
mean sky image, obtained from the actual data frames 
with sources masked out. Images were aligned by cen- 
troiding on objects in the field, and then combined with 
a clipped mean. 

Because there was significant cirrus during the run, 
zeropoints wer e set by measuring; the magnitudes from 
the FLAMEX (|Elston et alJl200ffl and WIRC ifs-band 
images. Because the Ks filter (central wavelength = 
2.15 jim, bandpass width = 0.31 urn) is only slightly nar- 
rower than K' (central wavelength = 2.12 jim, bandpass 
width = 0.35 fim), the magnitude differences between 
K' and Ks should be less than 0.1 magnitudes. As 0.1 
magnitudes is roughly our photometric uncertainty, we 
choose to report all measured magnitudes in the Ks fil- 
ter. 

3. MORPHOLOGY AND PHOTOMETRY 

Figure [4] shows the AO images of the three DOGs in our 
sample (left panels). Two of the three systems, DOGS 
1 and 2, are resolved. The third appears to be domi- 
nated by a point source, possibly the result of an AGN. 
Interestingly, there is no significant substructure in these 
images. The galaxy light profiles appear to be smooth, 
as opposed to comprised of multiple knots or nuclei. 

We used the 2D galaxy profile fitting routine GAL- 
FIT (Peng et al. 2002) to model each system. We pro- 
vided GALFIT with the image of each galaxy, its associ- 
ated PSF, and an initial guess for the best fitting galaxy 
model. GALFIT then matched for total magnitude, ef- 
fective radius, semi- minor to semi-major axis ratio, and 
position angle. We first ran eac h system with a simple, 
single Sersic model (jSersidll968h . 



TABLE 2 

Photometric and Morphological Properties of the 
DOGs 



E(r)=E c ca?p[-«((-) 1/n -l)], 



(1) 



allowing the Sersic parameter, L n\ and the effective ra- 
dius, r e , to float. The best fitting single Sersic parameter 
for each system is reported in Table [2j 

Images of the GALFIT derived models are shown in 
the middle column of Figure [H The right hand panels 
show the residuals after subtracting the model from the 
actual galaxy. DOG 1 is best fit by an exponential disk 
(actual best fit Sersic is n = 1.5) with an effective radius 
of 0.11". DOG 2 is best fit by a de Vaucouleurs profile 
(actual best fit Sersic is n = 3.5) with effective radius 
0.09". 

While both DOGs 1 and 2 are resolved they have very 
small effective radii. If they are at z = 1, 0.09" translates 
to a physical size of ~ 0.72 kpc. Because of the geom- 
etry of the universe, apparent size, for a given physical 
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2.92 


2.49 


1.15 


K s [mag] 


18.21 b 


18.42 b 


18.60 c 
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disk 


elliptical 


point source 
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a From Houck et al. 


(2005 


). b From 


FLAMEX 



(|Elstonet al] l2QQ6). c From WIRC. d A visual classifcia- 
tion from the AO images. 6 As measured by GALFIT f As 
measured by GALFIT g Concentration 11 Asymmetry 

size, remains roughly constant with redshift above z — \. 
That means that even if DOGs 1 and 2 are at z — 2, their 
physical sizes remain small. At z — 2, 0.09" translates 
to ~ 0.75 kpc. These are very small sizes, more typica l 
of today's blue compact dwarfs (e.g. JCorbin et aLl feOOG) 
than today's large star forming ULIRGs. 

DOG 3 is best fit by a point source of magnitude Ks ~ 
18.6. The structure seen in the image of DOG 3 appears 
to be primarily the result of structure in the AO PSF 
rather than "true" structure in the galaxy. This PSF 
structure is demonstrated in the GALFIT model of DOG 
3 (middle column) which is effectively a scaled version of 
the observed PSF (in this example, we used images of 
the tip-tilt star for DOG 2 as the PSF). 

After subtracting the best fitting model galaxies from 
the actual DOG images, the residual images (right hand 
column of Figure [4]) show little additional substructure; 
the possible exception being DOG 1. Any remaining flux 
within the residuals is small, with RMS deviations less 
than 5% of the ffux/pix at the centers of each DOG. DOG 
1 is the only object of the three suggestive of a secondary 
central component. We used GALFIT to simultaneously 
fit DOG 1 for both a disk and a bulge. The fit for the 
central component is actually disk-like with a Sersic in- 
dex of n = 0.3, similar to th e truncated central disk in 
NGC 2787 (|Peneret al.ll2002h . Although GALFIT pro- 
duced a fit for this second component the light within it 
was roughly a factor of 100 less than in the primary com- 
ponent, and the improvement over a single Sersic model 
was negligible. 

In order to compare the morphologies of the DOGs 
with other samples of high redshift dusty galaxies, we 
also measured two non-parametric morphological in- 
dexes, concentration (C), and asymm etry ( A). W e used 
the C and A definitions given in IConselicd (|2Q03h . C = 
5log (r so /r2o), where r$o is the radius that contains 80% 
of the light, and r2o is the radius containing 20% of the 
light. We measured the concentration index with circu- 
lar aperture photometry and a curve of growth technique. 
To provide a reference, we first measured the concentra- 
tion of the AO PSF, which has C ~ 5. While all three 
DOGs show high concentrations, 4 > C > 3, typical of 
early- type disks and e lliptical galaxies in the local uni- 
verse (|Conselicd [2003). they are less concentrated than 
the AO PSF. This is true even for the point-like DOG 3. 

Asymmetry is given by the normalized residuals after 
subtracting a 180-degree, rotated image of the galaxy 



5.0 
4.5 
4.0 



f 3.5 

II 3.0 
o 



2.5 

2.0 
0. 



i 



4 



.A A. 



"•- -#-- --©- - 
"• 



-O 
-A--0 



• # - 



■"G#" 



...::8 a 

-o o 



-O 0-9 



0.1 0.2 0.3 0.4 0.5 

Asymmetry 



Fig. 5. — Asymmetry (A) and concentration (C) measurements 
of dusty galaxies. Filled symbols are asymmetry measurements 
made in the X'-band. Open symbols are asymmetry measure- 
ments made in z-band. Measurements of the DOGs are shown as 
stars, and are labeled. A and C measurements for a sample of 
z ~ 2 SCUBA sources measured in HST ACS z-band images are 
shown as triangles. A second comparison sample of z ~ 1 LIRGs 
observed in both the K' and z-band are shown as circles. Lines are 
drawn between the two A measurements of each LIRG, with con- 
centration measures made in the X'-band. The C and A measures 
for the two AO PSF stars are also shown (squares). The spatially 
resolved DOGs show low asymmetry compared with the SCUBA 
sources. The LIRGs show similarly low asymmetry in the K' band 
images. Interestingly, the LIRGs show much higher asymmetry in 
the z-band (rest-frame blue-UV). The DOGs tend towards higher 
concentrations than the comparison samples. 

from itself. For a complete desc ri ption of the asymmetry 
parameter, see IConselice et al.l (J2000). Because asym- 
metry measures ar e dependent on identifying the galaxy 
center, we take the lConselicd (J2003) suggestion to use the 
minimum A after using a grid of rotation centers. Not 
surprisingly, the smooth images of DOGs 1 and 2 show 
small asymmetries, A < 0.1, typical of non-i nteracting 
galaxies in the local universe (IConseliceJl2003[ ). DOG 3, 
which is point-source dominated, has a higher A measure 
than the other two, A=0.39. Some portion of the asym- 
metry measured in DOG 3 can be attributed to struc- 
ture within the AO PSF. This structure can be seen in 
the GALFIT model image for DOG 3 (Figure gj). The 
asymmetry measures of the AO PSFs are A ~ 0.23. The 
non-parametric morphology measurements are summa- 
rized in Table El 

The uncertainties in these morphological measures are 
comprised of two distinct components, photometric un- 
certainties and PSF uncertainties. The high signal-to- 
noise ratio images of DOGs 1 and 2 translate into low 
photometric uncertainties. Because DOG 3 was observed 
for significantly less time, photometric uncertainties are 
more significant. To quantify the affect of photomet- 
ric uncertainties on our morphological measurements we 
have populated the AO images with model galaxies us- 
ing the IRAF ARTDATA package. These models were 
created to have similar morphological properties as the 
actual DOGs, but are placed at 10 different locations in 



the AO frames. For DOGs 1 and 2 morphological mea- 
sures of the 10 models show small dispersion, typically 
5% or less for both the GALFIT and C and A measures. 
For DOG 3 the dispersion reaches 20%. 

To quantify the contribution of PSF uncertainty to our 
overall error budget, we run the GALFIT routine with 
both the best guess PSF, and with a second PSF. For 
DOG 1, we use the PSF for DOG 2 as the secondary 
PSF. For DOG 2 we use the PSF for DOG 1. Since no 
PSF measurement was made for DOG 3 we use PSFs 1 
and 2 and compare. For DOGs 1 and 2, the change of 
PSF has a significant affect on the measured Galfit pa- 
rameters, changing them on the order of 20%. For DOG 
3 the change in the measured Galfit parameters from 
PSF 1 to PSF 2 is only 10%, less than the uncertainty 
from photometric noise. While we do not know the true 
uncertainties in the PSFs, these are good upper bounds 
on PSF affects on the Galfit morphological measures. 

4. DISCUSSION 

To place the three DOGs in our sample into context, 
we compare their morphological parameters with two ad- 
ditional samples of intermediate to high redshift, dusty 
galaxies. The firs t is a sample of of z = 2 SCUBA sources 
(sub-mm sources. [Pope et al. 2005) observed in the Great 
Ob servatories Origins Dee p Survey north field (GOODS- 
N; iGiavalisco et all I20Q4D . Although SCUBA sources 
have a different selection criteria, they are also very 
dusty. The sub-mm emission is thought to arise from cool 
dust in very strongly star-bursting galaxies. The second 
is a sample of z ~ 1 LIRGs from GOODS-S observed in 
the if '-band with Keck AO by the Center for Adaptiv e 
Optics Treasury Survey ( CATS: [Melbourne et ai]l2008[ ). 
These two samples where chosen because they contain 
dusty galaxies, and because they span the redshift range 
of the DOGs. 

Figure [5] is a plot of asymmetry vs. concentration for 
the DOGs in our sample (stars). Also plotted are the 
GOODS-N SCUBA sources (triangles) and the GOODS- 
S LIRGs (circles). The C and A measures for the two AO 
PSF stars are shown as well (squares). For the SCUBA 
sources, we only plot C and A measurements of galax- 
ies that were reporte d to have high signal-to-noise ra- 
tio (|Pope et al.l 12005). The A and C measurements for 
SCUBA sources were made in the reddest of the GOODS 
HST ACS bands, the z-band. These images are signif- 
icantly bl uer than o u r if '- band images of the DOGs. 
However, iPope et all (|2005l ) suggest that there is little 
change in the A and C measurements through to the 
near-IR .fi-band (based on a small subset for which NIC- 
MOS ii"-band is available). They propose that the high 
asymmetry in the SCUBA sample indicates significant 
recent merging. 

Also plotted on Figure[5]are the A an d C measurements 
of a set of z ~ 1 LIRGs from CATS ([Melbourne et all 
120081 ). For the LIRGs, both K'- band (filled circles), 
and z-band (open symbols) asymmetries are shown. The 
if'-band asymmetries of the LIRGs are low in compar- 
ison wit h the SCUBA s o urces . This was anticipated 
because iMelbourne et al.l (|2Q08h demonstrated that the 
morphologies of these LIRGs are primarily disk-like and 
unlikely to have had a recent major merger. Interest- 
ingly, if we measure the asymmetries of the LIRGs in 
bluer bands, for instance the z-band (rest-frame opti- 
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Fig. 6. — Simulated AO images of mergers. Images are \" 
on a side. The mergers are based on the AO image of DOG 1. 
Two copies of the DOG 1 image are added together after shifting 
spatially and scaling to different flux ratios. The left column shows 
the image of the merger. The middle column is the best fitting 
single Sersic model, as given by GALFIT. The right hand column 
gives the residual difference between the two. For comparison, 
the top row shows the image of DOG 1 with no merger added. 
Successive rows show mergers with declining flux ratios, 1/1, 1/3, 
1/5, and 1/10. The merger components are separated by 0.2" or 
1.6 kpc at z = 2. This is well below the seeing limited resolution 
of Keck and at or below the spatial resolution of HST NICMOS in 
the near-IR. At these separations we can clearly identify mergers 
with flux ratios of 1/5 and larger. 

cal to UV), we find that they scatter towards higher 
asymmetries, similar to the SCUBA sources. The LIRGs 
are more asymmetric in the bluer bands because of 
widespread star formation and dust obscuration, rather 
than merging. This suggests that interpretation of high 
asymmetry is difficult, especially when looking at rest- 
frame blue to UV images. 

The two resolved DOGs in our sample have, on aver- 
age, lower asymmetry than the SCUBA sources. Their 
asymmetry measures are similar to the if '-band asymme- 
tries of the LIRGs. This may be evidence against recent 
merging within the DOGs. In order to quantify this we 
generate model merger systems and analyze them with 
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Fig. 7. — The best fit single Sersic models of the DOGs (stars). 
Also shown are the comparison sample z ~ 1 LIRGs (circles). The 
left panel shows effective radius as a function of Sersic index. The 
right panel shows effective radius as a function of Ks magnitude. 
The two resolved DOGs show significantly smaller (roughly by a 
factor of four) effective radii compared with LIRGs of similar Ser- 
sic index or magnitude. If we translate these apparent sizes into 
physical sizes using the range of possible DOG redshifts z = 1 — 3, 
then the physical sizes of the DOGs are also about a factor of 
four smaller than the LIRGs for the full range of redshifts. This 
is because apparent sizes are relatively fixed for a given physical 
size above redshift 1. The DOGs are fundamentally different from 
z = l LIRGs. 

GALFIT. The models are based on the image of DOG 
1. We add two copies of the DOG 1 image together af- 
ter shifting spatially and scaling to different flux ratios. 
Figure [6] shows model mergers generated with an 0.2" 
separation. This separation is well below the resolution 
of seeing limited images and at or below the resolution 
of HST NICMOS in the near-IR. As with Figure HJ the 
left column shows the merger images. The middle col- 
umn shows the best fitting single Sersic model for each 
merger, provided by GALFIT. The right hand column 
gives the residual difference between the two. For com- 
parison, the top row shows the image of DOG 1 without 
any merger added, and a single Sersic model subtracted 
off. Successive rows show mergers with decreasing flux 
ratios, of 1/1, 1/3, 1/5, and 1/10. Visually, at these 
separations, 1/5 ratio mergers are easily detected in the 
residual images, but 1/10 ratio mergers are not obvious. 
When we cut the separation in half to 0.1", 1/5 ratio 
mergers become more difficult to identify, but 1/3 ratio 
mergers are still relatively easy to detect. Thus the AO 
images rule out 1/5 ratio mergers at 0.2" separations, 
and 1/3 ratio mergers at 0.1" separations. If the DOGs 
in our sample are at z = 2, 0.1" translates to a physical 
separation of 0.8 kpc. This is a strong constraint against 
ongoing major mergers as the trigger for the DOGs' IR 
energy. Note, we can not rule out late stage mergers in 
which the nuclei have already merged completely. 

While the LIRGs and the resolved DOGs in our sample 
have similar if'-band asymmetries, the DOGs have, on 
average, higher concentrations. The DOGs also tend to 
be more highly concentrated than the SCUBA sources. 
We will discuss possible implications for the high concen- 
tration below, but first we compare the physical sizes of 
the DOGs and the z = 1 LIRGs. 

In addition to measuring non-parametric morpholo- 



gies, we have measured 2D galaxy profiles for the DOGs. 
As we showed in the previous section, the three DOGs in 
our sample show a range of morphologies from disk-like 
(DOG 1), to elliptical-like (DOG 2), to point-like (DOG 
3). Figure [3 shows the best fit single Sersic parameters 
for the DOGs in our sample (stars). We make the same 
measurements on the AO data for the z ~ 1 LIRGs (cir- 
cles). Apparent sizes for the two resolved DOGs are a 
factor of 4 - 5 smaller than for LIRGs with similar Ser- 
sic indices or if -band magnitudes. We do not know the 
redshifts for the three DOGs in our sample, therefore we 
cannot directly compare physical radii. We do, however, 
know that DOGs are typically at z ~ 2 and span a range 
from z ~ 1 to 3. If the DOGs in our sample happen 
to be at the same redshift as the LIRGs (z ~ 1), then 
their physical radii are smaller than the LIRGs by about 
factor of four. Above z = 1, apparent size remains rel- 
atively constant with redshift for a fixed physical size. 
This means that even if the DOGs are at z = 3 they still 
have sizes about a factor of four smaller than the LIRGs 
of z = l. 

In order to test whether the small sizes observed for the 
DOGs are the result of an observational bias against mea- 
suring objects with extended radii, we embedded model 
galaxies into the DOG frames. These galaxies were de- 
signed to have the same if -band fluxes as the DOGs, 
but with radii similar to the LIRGs. While these model 
objects were lower surface brightness than the actual 
DOGs, we were able to accurately measure their sizes and 
luminosities with GALFIT to within better than 10%. 
This argues that the DOGs really do have small effective 
radii. It is a further indication that the DOGs are not 
the same type of objects as the z — 1 LIRGs. 

So what are the DOGs? The DOGs have very red rest- 
frame UV to IR luminosities, a magnitude redder than 
Arp 220. This suggests that, not only do they contain sig- 
nificant amounts of dust, but the dust may be ubiquitous, 
blocking the majority of UV emitting sources. While our 
sample is currently very small, it may turn out that in 
order to be selected as a DOG, a galaxy must have a 
high concentration and small physical size. Otherwise 
the likelihood of observing some UV light, unobscured 
by dust, increases. This hypothesis will be interesting to 
explore with larger samples. 

Another key piece of information for explaining the na- 
ture of the DOGs may be their SEDs. At least two of 
the DOGs show power-law SEDs suggestive of AGN ac- 
tivity. Oddly enough those two DOGs do not show a 
significant central point source in the AO image. Both 
of these sources, however, have more highly concentrated 
light profiles than the LIRGs. This may mean that they 
contain a dust enshrouded AGN, increasing their con- 
centration, even though the AGN is not obvious in the 
images. Based on a stacking analysis of v ery deep x- 
ray da ta in the Chandra Deep Field South iFiore et all 
(2008) show that DOGs in GOODS-S are likely to con- 
tain dust obscured AGN. DOG 3 does show point-like 
morphology suggesting that it too may contain a strong 
central AGN. Surprisingly, DOG 3 shows the least power- 
law like SED. The "bump" in the SED of DOG 3 may 
be indicating that there is also strong star formation in 
this galaxy, and while the morphology is point like, the 



C and A measures for DOG 3 suggest that it is not a 
pure point source galaxy. Clearly we need to increase 
the sample size in order to better understand the rela- 
tionship between morphology and SED type in the DOG 
population. 

In addition to morphology, spatially resolved spec- 
troscopy and kinematics may help place further con- 
straints on AGN activity or mergers in the DOG sample. 
It may be possible to obtain such data with AO-fed in- 
tegral fie ld unit (IFU) sp ectrographs on large telescopes 
(e.g.. IWright et al1l2007l ). Spatially resolved, rest-frame 
optical spectroscopy targeting Ha and [ N II] can re- 
veal t he presence and location of AGN (e.g. jBrand et all 
l2QQ7h . In addition, IFU spectroscopy can provide a mea- 
sure of the kinematics within a galaxy. While one of our 
systems has a disk-like light profile, we do not know if it 
is a normal rotating disk. Spatially resolved spectroscopy 
might reveal that this disk-like system is actually kine- 
matically disturbed in some way. 

We have demonstrated the power of high spatial res- 
olution AO imaging for measuring the morphologies of 
dust obscured galaxies beyond z — 1. We plan to con- 
tinue imaging the DOGs with Keck AO, and hope to add 
kinematic information with AO-fed integral field spec- 
trographs. These future observations may untangle the 
nature of the DOGs and help explain the formation of 
today's massive galaxies. 
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